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We report measurements of the magnetic penetration depth A in single crystals of CeCoIns down 
to ~0.14 K using a tunnel-diode based, self-inductive technique at 28 MHz. While the in-plane 
penetration depth tends to follow a power law, A// ~ y 3 / 2 , data are better described as a 
crossover between linear (T > T*) and quadratic (T <C T*) behavior, with T* the crossover 
temperature in the strong-coupling limit. The c-axis penetration depth Aj_ is linear in T, providing 
evidence that CeCoIns is a d-wave superconductor with line nodes along the c-axis. The different 
temperature dependences of A// and X± rule out impurity effects as the source of T*. 



The compounds CeMIns (M = Co, Ir, Rh) have re- 
cently been added to the heavy-fermion family, and have 
attracted much interest due to their similarity with the 
cuprates: quasi-2D structure and proximity to mag- 
netic order Qj. CeCoIn 5 , in particular, is a good can- 
didate for study: its superconductivity is not sensitive to 
small changes in unit-cell volume or composition, unlike 
CeCu2Si2, and it has the highest T c (~2.3 K) among the 
heavy-fermion superconductors. CeCoIns has tetragonal 
HoCoGas crystal structure, consisting of alternating lay- 
ers of Celri3 and 'C0I112' El . De Haas- van Alphen (dHvA) 
data revealed that the Fermi surface (FS) is quasi-2D, 
with an open 2D undulating cylinder extending along 
the [001] direction, as well as the large effective masses 
of electrons on this FS |!| . 

Recently, there has been mounting evidence for un- 
conventional superconductivity in CeMIns . Specific heat 
data reveal a T 2 term at low temperature, consistent with 
the presence of line nodes in the superconducting en- 
ergy gap ||. Thermal conductivity measurements with 
in-plane applied field show four-fold symmetry, consis- 
tent with nodes along the (±7r, ±7r) positions j|. NQR 
measurements show that there is no Hebel-Slichter peak 
just below T c H. Below T c the spin susceptibility is 
suppressed, indicating singlet pairing j|, |(|. However, 
there are some ambiguities in some of the measurements. 
Thermal conductivity data yield a X 837 low-temperature 
behavior, that the authors claim is close to T 3 behav- 
ior predicted for unconventional superconductors with 
line nodes in the clean limit ||. NQR measurements 
did not show the T 3 low-temperature behavior of 1/ Tj 
that is expected for a line node gap; instead 1/Tj sat- 
urates below 0.3 K |5f|. Microwave measurements down 



to ^0.2 K showed a non-exponential behavior, and the 
authors claimed that A(T) ~ T below 0.8 K @, though 
the data clearly show some curvature in that tempera- 
ture range. Further, the field was applied along the ab- 
plane, so the shielding currents have both in-plane and 
inter-plane components. In this paper, we present high- 
precision measurements of in-plane A// and inter-plane 
Aj_ penetration depths of CeCoIns at temperatures down 
to 0.14 K. We find that A// is best treated as a crossover 
from ~ T to ~ T 2 at a temperature T* . Combined with 
the result that A^ oc T, this gives strong evidence for 
non-local behavior in a d-wave superconductor as pre- 
dicted by Kostzin and Leggett |§. 

Details of sample growth and characterization are de- 
scribed in Refs. jj], [|. Measurements were performed 
utilizing a 28 MHz tunnel diode oscillator (nj with a 
noise level of 1 part in 10 9 and low drift. The magni- 
tude of the ac field was estimated to be less than 5 mOe. 
The cryostat was surrounded by a bilayer Mumetal shield 
that reduced the dc field to less than 1 mOe. The sample 
was aligned inside the probing coil in two directions: (1) 
ab plane perpendicular to the rf field, measuring the in- 
plane penetration depth X// (screening currents in the ab 
plane); or (2) with the rf field parallel to the plane, giving 
a combination of A// and Aj^. The sample was mounted, 
using a small amount of GE varnish, on a rod made of 
nine thin 99.999% Ag wires embedded in Stycast 1266 
epoxy. The other end of the rod was thermally connected 
to the mixing chamber of an Oxford Kelvinox 25 dilution 
refrigerator. The sample temperature is monitored using 
a calibrated Ru02 resistor at low temperatures (Tf, ase - 
1.8 K), and a calibrated Cernox thermometer at higher 
temperatures (1.3 K - 2.5 K). We report data only for 
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FIG. 1: Low-temperature dependence of the in-plane pene- 
tration depth AA//(T). Lower inset shows AA//(T) over the 
full temperature range. Upper inset shows AA//(T) vs T 1 ' 8 
in the temperature range (0.14-1.13) K. The solid line is a 
guide to the eye. 



T > 0.14 K. The value of T c was determined from mag- 
netization measurements to be 2.3 K, identical to the 
previously reported value [3]. 

The deviation AA(T) = X(T) - A(0.14 K) is pro- 
portional to the change in resonant frequency Af(T), 
with the proportionality factor G dependent on sample 
and coil geometries. For a square sample of side 2w, 
thickness 2d, demagnetization factor N, and volume V, 
G is known to vary as G oc i?3_o(l-iV)/ V, where R^d 
= w/[2(l+(l+2d/Jf)a,rctan(w/2d)-2d/w] is the effective 
sample dimension For our sample 2w k 0.73 mm and 
2d m 0.09 mm. We determine G from a single-crystal 
sample of pure Al by fitting the Al data to extreme non- 
local expressions and then adjusting for relative sample 
dimensions. Testing this approach on a single crystal 
of Pb, we found good agreement with conventional BCS 
expressions. 

Fig. shows AA//(T) as a function of temperature. 
We see that A A //(T) varies strongly at low temperatures, 
inconsistent with the exponential behavior expected for 
isotropic s-wave superconductors. On the other hand, 
the variation is not linear, but has an obvious upward 
curvature, unlike the low-temperature behavior expected 
for pure d-wave superconductors. A fit of the low tem- 
perature data to a variable power law, AA//(T) = a + 
bT l yields n = 1.43 ± 0.01 for sample 1 and 1.57 ± 0.01 
for sample 2. The upper inset of Fig. [j] shows this ap- 
proximate T 8 / 2 behavior for sample 1. Kosztin et. al. 
jl2], [l3[ [u| have proposed a theory that gives a T 3/2 
term from the gradual evolution of the pseudogap above 
T c to the superconducting gap below T c . While resistiv- 
ity measurements suggest the possibility of a pseudogap 
in CeCoIns jl5|], which renders this interpretation feasi- 
ble, a decrease in Knight shift was observed only starting 
at T c B. We take the latter to rule out a pseudogap 



mechanism. 

Before considering novel excitation processes, we note 
the important distinction between AA(T), which is di- 
rectly measured, and the superfluid density [p( T) — 
A 2 (0)/A 2 ( X)] which can be inferred only with the knowl- 
edge of A(0) |ll|. In the d-wave model, even if p 
varies strictly with T, i.e. p—l-aT/T c , the penetration 
depth is non-linear: X(T) = A(0) [1+1/2 (aT/T c ) + 3/8 
(a T/T c ) 2 + ...]. Hence there is always a quadratic com- 
ponent to A whose strength depends on a, which in the 
d-wave model, is inversely proportional to dA(9)/d9\ no d e , 
the angular slope of the energy gap at the nodes M. If 
p(T) is linear in T, there is no need to invoke another 
mechanism. 

To extract the in-plane superfluid density from our 
data, we need to know A// (0) . For a quasi-2D supercon- 
ductor with a cylindrical Fermi surface and the material 
parameters in Ref. || Q , we obtain A// (0) = 2600 A, con- 
siderably larger than the experimentally obtained value 
of 1900 A Q. This along with a large heat-capacity jump 
at T c leads us to consider strong-coupling corrections as 
listed below fl M: 



r^M = l + 1.8(^) 2 (ln(^) + 0.5); 



??A M = l + 5.3(^) 2 ln(^); 
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each r\ represents the correction to the corresponding 
BCS value. If we take the experimental value of AC/7 T c 
= 4.5 then Eq. [l] gives the characteristic (equiv- 
alent Einstein) frequency ojo = 9.1 K and A^y(0) = 

1500 A. However, Petrovic et. al. JlJ argued that since 
C/T increases with decreasing temperature, the specific 
heat coefficient 7 is temperature-dependent below T c . 
This effect calls into question simple estimates of strong- 
coupling corrections for CeCoIn 5 . A better estimate is 
to use AC /AS, where AS is the measured change in 
entropy of the sample from T = to T c . Ref. [I] then 
gives AC/ AS — 2.5, so that ojo — 17.9 K, resulting in 
Ag c = 2.1k B T c and A^(0) = 2000 A. On the other hand, 
the larger AC of Ref. || yields AC / AS = 2.8 and uj 
= 15 A K, leading to Ag c = 2.2k B T c and A^(0) = 1900 

A. These values of Ayy(0) are close to that obtained by 
Ormeno et. al. 0. 

Although we will argue that non-local effects are im- 
portant, we will refer to (A//(0)/A//(T)) 2 as the "su- 
perfluid density." Fig. [2] shows the calculated behavior 
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FIG. 2: Low-temperature in-plane superfluid density pt/(T) 
= [A//(0)/A//(T)] calculated from AX /f (T) data in Fig. 1 
(thick lines). The thin lines correspond to fits to data using 
Eq. H using three values of A//(0). Inset shows p//(T) over 
the full temperature range. 
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FIG. 3: Low-temperature dependence of inter-plane (open 
circles) penetration depth AX±(T), after subtracting the in- 
plane component. In-plane AA//(T) (crosses) data is also 
shown for comparison. Solid line is a linear fit from 0.14 K to 
1 K. Inset shows inter-plane superfluid density p±(T) for the 
whole temperature range. 



of that quantity using the three values of X// (0) ob- 
tained above. We follow the procedure in Ref. [U] to 
compute the experimental superfluid density, using the 
T 8 / 2 fit to estimate the small difference between A//(0) 
and X// (0.14 K). In each case, p{T) is clearly not linear 
in T. 

Non-linearity in p( T) can arise from a crossover from 
an intermediate-temperature (pure) linear- T behavior 
to, for example, low-temperature (impurity-dominated) 
quadratic behavior as pointed out by Hirschfeld and 
Goldenfeld (23]. They interpolated between these two 
regions using 

X = X Q + bT 2 /(T*+T), (4) 

where T* is the crossover temperature. In terms of su- 
perfluid density, one obtains |ll| 

. > aT 2 /T c 
A P//( T ) = J^pf> (5) 

where T* depends on impurity concentration. 

A much more provocative source of the crossover of 
Eq. H was suggested by Kosztin and Leggett [KL] ||, 
who showed that for rf-wave superconductors, nonlocal 
effects change the linear behavior to quadratic below a 
crossover temperature T* local = A o ^//(0)/A // (0). 

The solid lines in Fig. g are fits to Eq. || and are very 
good for all three values of A//(0). The value of a varies 
from —0.5 to 0.7, the smallest value of a belonging to the 
largest value of A//(0). The value of a obtained is sim- 
ilar to that found for YBa 2 Cu 3 6 .95 (a - 0.6) ||, pi, 
but smaller than that of Tl 2 Ba 2 Cu06+,5 (a ~ 1.0) ]29] 
and K-(ET) 2 Cu[N(CN) 2 ]Br (a - 1.2) [jf§. The value of 



T* varies less across the three A//(0) values, from 0.32 K 
to 0.42 K. These values of T/T c (-0.14 - 0.18) differ 
from the cuprates (23| H| |6) and the organic supercon- 
ductor K-(ET) 2 Cu[N(CN) 2 ]Br (-0.05), where impurity 
scattering is presumed to be the source. Further, Ref. ^ 
puts an upper limit of 20 ppm on the impurity concen- 
tration. In the dirty c?-wave model J22j, this gives the 
unitary-limit scattering rate T — 1.5 x 10 8 s _1 , which 
yields an upper limit for T* — 65 mK. This is about 
5 times smaller than the experimentally obtained values 
above, suggesting that the sample is too clean for the 
dirty d-wave model to be applicable. 

Having ruled out impurity scattering, we turn to non- 
local electrodynamics as the source of the crossover in 
pi / (T) . For a d-wave superconductor with line nodes 
along the c-axis, nonlocality is expected to be relevant 
only when the applied magnetic field is oriented parallel 
to the c-axis, while the effect of impurities should not de- 
pend on the orientation of the field. As KL noted, if T* 
is noticeably smaller for ff 1 c than for H//c we may 
conclude that the observed effect is due mainly to non- 
local electrodynamics and not to impurities. For H _L c, 
screening currents flow both parallel and perpendicular 
to the c-axis, mixing X// and A^ with the frequency 

shift given by^ = ^(^ + ^)0, where V is 
the effective coil volume and /o the resonant frequency 
with the sample absent. In order to extract Aj_ we sub- 
tract out the Xj i component from A/j_. Fig. || shows 
the inter-plane penetration depth A_l of CeCoIns down 
to 0.14 K. It is clearly linear in T from 0.14 K to 1 K. 
To obtain the superfluid density, we estimate Aj_(0) from 
the H C 2 ani sotropy of —2.3 pi] ], and the fact that A(0) 
oc yjH c2 (0) 0' obtaining A ± (0) - 2700 A. This is close 
to the value of —2800 A obtained from microwave mea- 
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surements in the planar geometry p8[| . If we fit X±(T) 
to Eq. |[ we find T*^ < 0.15 K, significantly smaller than 
0.32 K obtained for the in-plane case. This satisfies the 
Kosztin-Leggett test and indicates that the superfluid 
response of CeCoIns is governed by nonlocal electrody- 
namics. This is also strong evidence that CeCoIn$ is a 
d-wave superconductor with line nodes along the c-axis. 
Sr 2 Ru04 failed this test |n| because its line nodes are 
horizontal instead of vertical. Kusunose and Sigrist ar- 
gued that horizontal line nodes give power-law behaviors 
with less angular dependence for any inplane direction of 
the screening currents, and hence applied field (^9|. A 
calculation of p± is shown in the inset of Fig. 0: the up- 
turn below 0.5 K is an artifact of the choice of A c (0). A 
larger value of A c (0) would remove this feature, but there 
is no justification for doing so. 

As a final test of the non-local scenario, we esti- 
mate Ty J using strong-coupling parameters. From the 
measured Jf C 2(0)[001] value of 49.5 kOe, the coherence 
length £//(0) is calculated to be 82 A Q. Together 
with the earlier-derived values of Ag c = 2.2A;b T c and 
A^(0) = 1900 A, we find the strong-coupling nonlocal 
crossover temperature T* nonlocal = Ag c £//(0)/A^(0) = 
0.22 K. Using a weak-coupling d-wave A(0) = 2.14fc j BT c , 
we find T^ onlocal = 0.26 K. We regard either value to be 
satisfactorily close to the experimental value of 0.32 K. 
Note that the value of £//(0) is different from the calcu- 
lated BCS value of 58 A_i or the strong-coupling cor- 
rected value of ~ 50 A jH This is not surprising since 
the BCS expressions fl9fl assume a spherical FS, while 
LDA band structure reveals a very complicated FS with 
contributions from three different bands |jo| . 

In conclusion, we report measurements of the magnetic 
penetration depth A in single crystals of CeCoIn 5 down 
to ^0.14 K using a tunnel-diode based, self-inductive 
technique at 28 MHz. The in-plane penetration depth 
(A//) exhibits a crossover between linear (at high T) and 
quadratic (low T) behavior with a crossover temperature 
^nonlocal ~ 0-^2 K. Such behavior can arise in a super- 
conductor with nodes in the gap either in a dirty ci-wave 
model or from nonlocal electrodynamics. The linear low- 
temperature dependence of the c-axis penetration depth 
Aj_ strongly favors the nonlocal model with line nodes 
parallel to the c-axis. We also demonstrate that strong- 
coupling corrections are required to reconcile various ex- 
perimentally determined superconducting parameters. 
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